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We study the coupling parameter dependence of the Higgs boson mass in a chirally invariant 
lattice Higgs- Yukawa model emulating the same Yukawa coupling structure as in the Higgs-fermion 
sector of the Standard Model. Eventually, the aim is to establish non-perturbative upper and lower 
Higgs boson mass bounds derived from first principles, in particular not relying on vacuum stability 
' considerations for the latter case. Here, we present our lattice results for the lower Higgs boson mass 

bound at several values of the cutoff A and compare them to corresponding analytical calculations 
Mh , based on the effective potential as obtained from lattice perturbation theory. Furthermore, we give 

' a brief outlook towards the calculation of the upper Higgs boson mass bound. 

(N ■ 



I. INTRODUCTION 



With the existing evidence for the triviality of the Higgs sector P, H, H, 0, H, H, 0] of the Standard 
Model, rendering the removal of the cutoff A from the theory impossible, physical quantities in this 
^ ' sector will, in general, depend on the cutoff. Though this restriction strongly limits the predictive 

power of any calculation performed in the Higgs sector, it opens up the possibility of drawing conclu- 
sions on the energy scale A at which new physics has to set in, once, for example, the Higgs boson 
mass has been determined experimentally. 

Besides the obvious interest in narrowing the energy interval of possible Higgs boson masses consis- 
I tent with phenomenology, the latter observation was the main motivation for the great efforts spent 

on the determination of cutoff-dependent upper and lower Higgs boson mass bounds. In perturbation 
Tlj- ' theory such bounds have been derived from the criterion of the Landau pole bein g si tuated beyond 

the cutoff of the theory (see e.g. [sl. [ol. [lOj ) . from unitarity requirements (see e.g. flll . lT2l .fi3j) and from 
vacuum stability considerations (see e.g. [1, fH, fH, fH, fl3j fi3| ) , as reviewed in |l9f . 
^\ \ However, the arguments relying on the vacuum instability have to be considered with care. It is 

■ well known that the constraint effective potential, which becomes equal to the effective potential in 

the thermodynamic limit, is convex and can therefore never be unstable |i2Q] . Additionally, it has 
recently been argued that the appearance of vacuum instability is actually caused by the incorrect 
removal of the cutoff from the trivial theory [U, fg^]. In this work, we will not address the question 
of the validity of the vacuum instability argument but rather investigate whether the lower mass 
constraint can also be established without relying on vacuum stability considerations. Additionally, 
the validity of the perturbatively obtained upper Higgs boson mass bounds is unclear, because the 
corresponding perturbative calculations had to be performed at rather large values of the renormalized 
quartic coupling constant. The latter two remarks make the Higgs boson mass bound determination 
an interesting subject for non-perturbative investigations starting from first principles, such as the 
lattice approach. 

The main objective of lattice studies of the pure Higgs and Higgs- Yukawa sector of the electroweak 
Standard Model has therefore been the non-perturbative determination of the cutoff-dependence of the 
upper and lower Higgs boson mass bounds 0, fE flj flllj fl2, fill as well as its decay properties [24|. There 
are two main developments that warrant the reconsideration of these questions. First, with the advent 
of the LHC, we are to expect that properties of the Standard Model Higgs boson, such as the mass 
and the decay width, will be revealed experimentally. Second, there is, in contrast to the situation 
of earlier investigations of lattice Higgs- Yukawa models 0, fH, fS^, fH, fH, fM, fHI, which suffered 
from their inability to restore chiral symmetry in the continuum limit while lifting the unwanted 
fermion doublers at the same time, a consistent formulation of a Higgs- Yukawa model with an exact 
lattice chiral symmetry [s^l based on the Ginsparg- Wilson relation fsSj . This new development allows 
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to maintain the chiral character of the Higgs-fermion coupHng structure of the Standard Model on 
the lattice while simultaneously lifting the fcrmion doublcrs, thus eliminating manifestly the main 
objection to the earlier investigations. The interest in lattice Higgs- Yukawa models has therefore 
recently been renewed [s^, IH, [H, |37| . 

Before, however, questions of the Higgs boson mass bounds and decay properties can be addressed, 
the phase structure of the model needs to be investigated in order to determine the (bare) coupling con- 
stants in parameter space where eventual calculations of phenomenological interest can be performed. 
Much effort has therefore bee n sp ent on phase structure investigations of lattice Higgs- Yukawa models 
in the past, see e.g. Refs. [3l,[3l^[43,lll|,|42,[4l,|4l,|4^|4^ for an only partial hst. The phase structure 
of the new, chiraUy invariant Higgs- Yukawa model has been studied analytically by means of a large 
Nf computation (3^. [3^. where Nf denotes the number of fermion generations, as well as numerically 
by direct Monte-Carlo simulations [sslls^. 

In the present paper we study the dependence of the Higgs boson mass on the model parameters 
by direct Monte-Carlo simulations. We confirm that the smallest and largest Higgs boson masses are 
obtained at vanishing bare quartic self-coupling and at infinite bare quartic coupling, respectively, 
as expected from perturbation theory. We then present our results on the cutoff-dependence of the 
lower Higgs boson mass bound and examine finite volume effects. The resulting finite volume Higgs 
boson mass bounds are extrapolated to the infinite volume limit. Since the aforementioned results 
were determined in the mass degenerate case, i.e. with equal top and bottom quark masses, we also 
investigate the effect of the top-bottom mass splitting on the Higgs boson mass. All these numerical 
findings are compared to corresponding analytical calculations based on the effective potential at the 
1-loop order of lattice perturbation theory. The very good observed agreement between these two 
approaches justifies deriving the lower mass bound analytically from the effective potential also in the 
numerically very demanding, physical situation, where the bottom quark is approximately 40 times 
lighter than the top quark. Finally, we give a brief outlook towards the upper Higgs boson mass 
bound. 



II. THE SU(2)l X \J{1)r LATTICE HIGGS- YUKAWA MODEL 

The model we consider here, is a four-dimensional, chirally invariant SU{2)l x U{1)r lattice Higgs- 
Yukawa model [3^ . aiming at the implementation of the chiral Higgs-fermion coupling structure of 
the pure Higgs- Yukawa sector of the Standard Model reading 

Ly = Vb [t, b)^^bR + yt {t, b)^if>tii + c.c, (1) 

with yt^b denoting the bare top and bottom Yukawa coupling constants. Here, we have restricted 
ourselves to the consideration of the top-bottom doublet {t,b) interacting with the complex Higgs 
doublet If [ip = iT2if*, and are the Pauli-matrices), since the Higgs dynamics is dominated by the 
coupling to the heaviest fermions. For the same reason we also neglect all gauge fields in this approach. 

The fields considered in this model are one four-component, real Higgs field <i>, being equivalent 
to the complex doublet (p of the Standard Model, and Nf top-bottom doublets represented by eight- 
component spinors t/j^^^ = (t*^*-' , 6^'^), i — l,...,Nf. 

The chiral character of the targeted coupling structure in Eq. ([T|) can be preserved on the lattice 
by constructing the fermionic action Sp from the (doublet) Neuberger overlap operator V^"""^ 
according to 

Sf = (2) 

i=l 

M = + P+0t diag {yt , y,) P+ + P_ diag {yt , yb) ^P- . (3) 

Here the Higgs field $3,, defined on the site indices x — {x, t) of the x Lj-lattice, was written as a 
quaternionic, 2x2 matrix (f>x = $^cr^, ctq = 1, cr^ = —iTj with r denoting the vector of Pauli matrices, 
acting on the fiavour index of the fermionic doublets. The left- and right-handed projection operators 
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P± and the modified projectors P± are given as 



with p being the radius of the circle of eigenvalues in the complex plane of the free Neuberger overlap 
operator 47] . This action now obeys an exact SU(2) ^ x 11(1)^ lattice chiral symmetry. For e SU(2) 
and Un G U{1) the action is invariant under the transformation 

^ ^ UrP+^j + nLP-4', ^p^^pp+n[ + i,p^ ul, (5) 

(j^^URcj^nl, cl)^ ^^Lci^^ul (6) 

Note that in the mass-degenerate case, i.e. jjt — jjbi this symmetry is extended to SU(2)i x SU(2)/;. In 
the continuum limit the symmetry given in Eq. ^ and Eq. ^ recovers the continuum SU(2)i x U/j(l) 
chiral symmetry and the lattice Higgs- Yukawa coupling becomes equivalent to the coupling structure 
in Eq. ([1]) when identifying 

'Px = C- ( , (p^= iT2ip% = C • ^2 , -^1 , and yt.b = ^ (7) 



for some real, non-zero constant C. Note also that in absence of gauge fields the Neuberger Dirac opera- 
tor can be trivially constructed in momentum space, since its eigenvectors and eigenvalues i''^{p), p G V 
are explicitly known. Here, V is the set of all lattice momenta and e = ±1 is the sign of the imaginary 
part of the eigenvalues coming in complex conjugate pairs. In this notation they are given as ^3^] 

V {p) = P + P-— / p =2^sm(p^), p =42^sm ( — 1 (8) 

where r denotes the coefficient of the Wilson term in the underlying Wilson Dirac operator. 
Finally, the lattice Higgs action 5$ is given by the lattice ^''-action 

S^^-kY.'^l [<^,+^ + <f>,_^] + ^ $t <i,^ + A ^ ($t $^ _ Nf) ' , (9) 



which is equivalent to the lattice action in continuum notation 



(10) 



with the bare mass m, the bare quartic self-coupling constant A, and V^^ denoting the lattice forward 
derivative in direction /i. The connection is established through a rescaling of the Higgs field and the 
involved coupling constants according to 

^, = V2.(^^o_,$3 j, A = ^, m= ^ , = (11) 

In the given model the expectation value ((/j) would always be identical to zero due to the symmetries 
given in Eq. ([5]) and Eq. (I6|). To study the mechanism of spontaneous symmetry breaking, one 
usually introduces an external current, which has to be removed after taking the thermodynamic limit, 
leading then to the existence of symmetric and broken phases with respect to the order parameter 
as desired. An alternative approach, which was shown to be equivalent in the thermodynamic 
limit [i^, H^llOl, is to use the global symmetries of the model to rotate each field configuration {}Px) 
after its generation in the Monte-Carlo process according to 

= C^^. (12) 
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with the SU(2) matrix U selected for each field configuration {(px) such that 





E,„rot 



(13) 



Here, we use this second approach. According to the notation used in Eq. pop . which already in- 
corporates a factor 1/2 in contrast to the standard notation, the relation between the bare vacuum 
expectation value of the Higgs mode v and the expectation value of (^''°* is then given as 



(14) 



In this approach the unrenormalized Higgs field hx and the Goldstone fields f/xjfxiS? '^^'^ directly be 
read out of the rotated field (/5^°* according to 



9x H 
v + h, 



(15) 



III. SIMULATION STRATEGY AND OBSERVABLES 



The general method we apply to determine the lower and upper Higgs boson mass bounds is the 
numerical evaluation of the whole range of Higgs boson masses that are attainable within our model 
at a fixed value of the cutoff in consistency with phenomenology. The latter requirement restricts the 
freedom in the choice of the bare parameters k, yt f,, A due to the phenomenological knowledge of the 
renormalized vacuum expectation value of the Higgs mode (vev), i.e. Vr/a = 246 GeV, and of the 
top and bottom quark masses, i.e. mt/a ~ 175 GeV and miy/a « 4.2 GeV, respectively. Here, mt, 
rrib denote the lattice masses and a is the lattice spacing. For our numerical simulations we use the 
tree-level relation 



yt,b 



(16) 



as a first guess to set the bare Yukawa coupling constants yt and j/b. The actually resulting quark 
masses are discussed in section |Vl Furthermore, the model has to be evaluated in the broken phase, 
i.e. at non- vanishing vacuum expectation value of the Higgs mode, v ^ Q, however close to the 
phase transition to the symmetric phase. We also use the phenomenologically known value of the 
renormalized vev to determine the lattice spacing a and thus the physical cutoff A according to 



246 GcV = — = 
a 



A = 



Gg{p 



P 



''Gp 



Zg 



(17) 



Here, the Goldstone propagator mass nicp is zero except for finite volume contributions and the 
renormalization constant is obtained from the Goldstone propagator Gg(S>^) in momentum space 
with fp' denoting the squared lattice momenta. The Higgs propagator mass niHp and the corresponding 
renormalization factor Zu are obtained in the same way through the relation 



Zh 



where the Goldstone and Higgs propagators in position space are given as 



Gg{x,v) 



1 



3 

■^^{9x9y 



and GH{x,y) = {hxhy 



(18) 



(19) 



The aforementioned Higgs propagator mass is connected via perturbation theory to the actual 
physical Higgs boson mass mn given by the corresponding energy eigenvalue of the Hamiltonian. It 
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can be obtained from the exponential decay of the Higgs time shce correlation function 



1 



6 • '^a.t+At) with Cff (Lt/2 > > 1) oc cosh^m// • {Lt/2 - At)J . (20) 

t=l x,y 



Correspondingly, we derive the physical top and bottom quark masses mt , mi, from the fermionic time 
slice correlators 



with f = t,b selecting the quark flavour. We remark here that the full all-to-all correlator as defined 
in Eq. ()21|) can be trivially computed by using sources which have a constant value of one on a whole 
time slice for a selected spinor index and a value of zero everywhere else. This all-to-all correlator 
yields very clean signals for the top and bottom quark mass determination. 

For a given cutoff A the aforementioned requirements of reproducing the vev and the quark masses 
still leave open a one-dimensional freedom in the model parameters, which can be parametrized in 
terms of the quartic self-coupling constant A. However, aiming at lower and upper Higgs boson 
mass bounds, this remaining freedom can be fixed, since it is expected from perturbation theory 
that the lightest Higgs boson masses are obtained at vanishing self-coupling constant A = 0, and 
the heaviest masses at infinite coupling constant A = cxd, respectively, according to the qualitative 
one-loop perturbation theory result for the Higgs boson mass shift 



Subsequently, it will be explicitly checked by direct lattice calculations that the lowest and highest 
Higgs boson masses are indeed obtained at vanishing and infinite bare quartic coupling constants, 
respectively. The lower Higgs boson mass bound will then be searched for in the weak self-coupling 
region, i.e. at A <C 1, while the upper bound will be determined from simulations performed in the 
strong coupling region, i.e. at A 3> 1. 

For the numerical evaluation of the model we have implemented a PHMC-algorithm [5l|, [H, [53| , 
allowing access to the physical situation of odd Nf. We make heavy use of the fact that no gauge fields 
are considered here. In this setup the Dirac operator is diagonal in momentum space and its eigenvalues 
and eigenvectors are explicitly known. In our approach we use a fast Fourier transformation [s^ l to 
switch between momentum space, where can be trivially applied, and position space, where we 

perform the multiplication with the field Besides that, it was crucial to reduce the auto-correlation 
times of the measured observables, which we achieved with the help of Fourier acceleration [ssl [56l|. 
We also remark here that the condition number of the fermionic matrix could greatly be reduced by 
preconditioning. However, the details of the algorithm will be discussed elsewhere. 

All results in the following have been obtained at Nf = 1 with degenerate Yukawa coupling con- 
stants, i.e. yt = Ub (unless otherwise stated), owing to the much higher numerical requirements at the 
realistic ratio yb/vt = 0.024 exceeding our available resources. However, on some small lattice sizes 
we could also investigate the dependence of the obtained Higgs boson masses on the top-bottom mass 
splitting ratio yb/yt, allowing ultimately to extrapolate our results to the physical scenario. Addition- 
ally, we are currently working on the Nf = 3 results to account for the colour index (even though 
gauge fields are stfll absent). 

IV. A-DEPENDENCE OF THE HIGGS BOSON MASS AND EFFECTIVE POTENTIAL 

We now turn to the determination of the cutoff dependent lower Higgs boson mass bound. As 
discussed in the previous section, one would expect the lightest Higgs boson masses at vanishing 
quartic self-coupling parameter A = 0. However, one should remark here that the given argument 
is not complete, since the phase transition line changes while varying A as well, thus making the 
bare Higgs boson mass m a function of A for fixed cutoff A and fixed Yukawa coupling constants. In 
fact, the phase transition line is strongly shifted when varying the self-coupling parameter [33 |. From 




(21) 




(22) 
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Eq. (|22p alone one can therefore not conclude whether the lightest Higgs boson mass in this model 
will actually be observed at A = or rather at slightly stronger self-coupling. 

To examine this we consider the effective potential ^571] in terms of the vacuum expectation value v. 
In the large iV/-limit with yt.b ~ 0{1/ ^jNf) and A = 0{1/Nf) it was given for this model in Ref. [13] 

as 



U[v] 
UfM 



1 



m^v'' + Xv'^ + Uf[v] 



^2N 



f 



v+{p) +Vtv[l 



2p 



log 



+ybv(l 



2p 



(23) 



(24) 



For a given value of v, and thus for a fixed cutoff A, the dependence of the bare Higgs boson mass on 
the quartic self-coupling constant can be trivially determined by setting the first derivative of U[v] to 
zero according to 



m 



-4Au^ 



1 d 

V dv 



Uf[v] 



(25) 



Correspondingly, one can also estimate the Higgs propagator mass from the curvature of the effective 
potential yielding 



''Hp 



^8Xv^~-^Uf[v] + ^Uf[v] 



V dv 



dv'^ 



(26) 



From this result it is clear that the lightest Higgs boson masses should be observed at vanishing quartic 
coupling constant, since the fermionic contribution Uf does not depend on A at the considered order. 
This observation allows to restrict the search for the lower Higgs boson mass bounds to the setting 
A = in the following. 

We check the latter result by studying the A-dcpcndence of the Higgs boson mass in direct Monte- 
Carlo simulations performed on a 16^ x 32-latticc for constant cutoff A — 400 GeV, Nf = 1, and 
degenerate Yukawa coupling parameters fixed according to Eq. (fT6| . In Fig. [1^ we show the Higgs 
correlator mass m^r/a, which is almost identical to the propagator mass in the weak coupling region, 
versus the quartic coupling constant A. The mass monotonously falls with A going to zero and 
converges to the Higgs boson mass obtained at A = 0. We remark that at A = the model is still 
well defined even in the broken phase, at least sufhciently close to the phase transition line thanks to 
the squared bare mass being positive provided that the Yukawa coupling constants are non-zero. 
Furthermore, the numerical results are in good agreement with the analytical predictions in Eq. (|26p . 
where the renormalization constant Zq has been set to one for the conversion between A and v. 

In Fig. [1)3 we additionally present the numerical results for the mass shifts Sm'jj together with the 
corresponding squared bare masses . One sees that the bare mass becomes smaller with increasing 
self-coupling constant A, but that this effect is overcompensated by the much stronger rise of the mass 
shift Smfj, as expected from the predictions in Eq. (|25p and Eq. ([^5]) marked by the dotted lines. 
However, on a quantitative level there is a huge discrepancy between the analytical predictions and 
the numerical results concerning both, the mass shifts as well as the bare masses. 

This discrepancy can be fixed by taking the next order in the loop expansion of the effective 
potential [s^l of the lattice Higgs action in Eq. PH)) into account, which would correspond to the 
l/A^/-correction of the pure bosonic part of the effective potential in the language of the large Nf 
analysis. This next to leading order contribution is given by connecting two legs of the four-point 
vertex and summing up the arising scalar propagator over all lattice momenta yielding the improved 
effective potential 



imp 



[v] = U[v]+ Xv^ 



1 



6 



E 



6 



p^V 



r+m^p O^^iT'^ +"'Gp 



(27) 



Here, the momentum p = is excluded in the sum over the Goldstone loop, since the Goldstone 
propagator at zero momentum is identical to zero due to the rotation of the field introduced in 
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FIG. 1: (a) The dependence of the Higgs boson mass mn/a on the quartic self-coupling constant A at 
A = 400 GeV on a 16^ x 32-lattice for constant, degenerate Yukawa coupling constants fixed according to 
Eq. (|16|> and Nf = 1. The prediction from the effective potential in Eq. ((23} is represented by the dashed 
curve, (b) The corresponding Higgs boson mass shifts 5m?fj and bare masses versus the quartic coupling 
constant A. The dotted (dashed) lines show the analytical results derived from the unimproved (improved) 
effective potential. 



section ini However, with increasing lattice volume this exclusion looses any significance. The factors 
appearing in the nominators of Eq. I|27p are given by the multiplicities of the considered diagrams. 
Though the effective potential U[v\ was actually given in bare constants, we plug the renormalized 
Higgs and Goldstone boson masses into the expressions for the scalar propagators, in order to further 
improve the result. For the evaluation of Eq. ([?f|) independent of any knowledge from Monte-Carlo 
simulations we use togp = in the following and determine mup self-consistently with respect to the 
curvature of the resulting improved effective potential at its minimum. 

The results for the mass shifts and bare masses obtained by the improved effective potential are 
shown in Fig. [T}d depicted by the dashed lines. The improved results are in very good agreement 
with the numerical data. The reason why the unimproved effective potential could describe the 
renormalized masses presented in Fig. [1^ so well - in contrast to the bare masses and mass shifts - is, 
that the improvement term exactly cancels in Eq. (|26p but not in the equation for the bare mass. 



V. LOWER MASS BOUNDS 



Given the knowledge about the A-dependence of the renormalized Higgs boson mass we can now 
safely determine the cutoff-dependence of the lower Higgs boson mass bound by evaluating the Higgs 
boson mass at A = for several values of A. Two restrictions limit the range of accessible energy 
scales: on the one side all particle masses have to be small compared to A to avoid unacceptably large 
cutoff-effects, on the other side all masses have to be large compared to the inverse lattice size to avoid 
finite volume effects. As a minimal requirement we demand here that all particle masses m in lattice 
units fulfill 

m < 0.5 and rh ■ Ls,t > 2. (28) 

For a lattice with side lengths Lg — Lt — 32, a, degenerate top/bottom quark mass of 175 GeV, and 
Higgs boson masses ranging from 40 to 70 GeV one can thus access energy scales A from 350 GeV to 
approximately 1100 GeV. 

In Fig. [2^ we show the numerically obtained Higgs correlator masses mu/a versus the cutoff A. 
All presented results have been obtained in Monte-Carlo simulations with identical, degenerate bare 
Yukawa coupling constants fixed according to Eq. p^ . A = 0, and Nf = 1. To illustrate the influence 
of the finite lattice volume we have rerun all simulations with exactly the same parameter settings but 
different lattice sizes. Those results belonging to the same parameter sets are connected by dotted lines 
to guide the eye. Additionally, we also present the corresponding analytical, finite volume expectations 
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for the cutofF-dependence of the Higgs boson mass as derived from the effective potential. We observe 
that the effective potential calculation is in very good agreement with the cutoff dependence of the 
measured Higgs boson masses. It also correctly predicts the observed finite volume effects. 

The cutoff dependence of the corresponding top quark masses is presented in Fig. [2)3. Those results 
belonging to the same bare parameters are connected by dotted lines for a better overview. The 
renormalized top quark mass, and thus the renormalized Yukawa coupling constant, decreases slightly 
with growing cutoff as expected in a trivial theory. In order to hold the fermion masses constant as 
desired one would have to fine-tune the bare Yukawa coupling parameters, giving rise to an additional 
but small contribution to the cutoff dependence of the Higgs boson mass. 

For an efficient tuning of the bare Yukawa coupling constants, however, a thorough analytical 
understanding of the observed behaviour of the fermion mass would be indispensable. We therefore 
compare the presented numerical findings for the top quark mass to the corresponding 1-loop lattice 
perturbation theory predictions, depicted by the dashed lines in Fig. [2)3 and observe good agreement, 
allowing one ultimately to hold the fermion masses constant in follow-up Monte-Carlo calculations. 
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FIG. 2: (a) The lower Higgs boson mass bound m'^'^ /a versus the cutoff A. (b) The cutoff dependence of the 
top quark mass rrit/a. To illustrate finite volume effects, simulations have been rerun with identical parameter 
sets but different lattice sizes. Runs with same parameter sets are connected via dotted lines to guide the 
eye. The analytical predictions for the respective lattice sizes derived from the effective potential are depicted 
by the dashed curves. All presented results have been obtained in Monte-Carlo simulations with identical, 
degenerate bare Yukawa coupling constants fixed according to Eq. (fT6)l . A = 0, and Nf = 1. 



The presented perturbative results are given by the tree-level fermion propagator plus the 1- 
Goldstone-loop and 1-Higgs-loop corrections according to 



hh{p) 

r(p) 



D{p)-^{p)\ , D{p)=V^''-\p)+vBoT{p) 



B^V(k)b-\k)B^V{p) 



1 



''Hp 



^ ^ B,r{k)D-\k)B,r{p) ^.^^^ t = p^k 



(29) 
(30) 



(31) 



Here is the squared lattice momentum of the relative momentum t = p—k as defined in Eq. ([8]) and 
(p) denotes the 8 x 8 doublet Dirac matr ix in momentum space for the given lattice momentum 
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p G V. We remark that the contribution from the tadpole diagram to the fermion propagator is 
exactly canceled by the renormalization of the vev induced by the tadpole diagram. Again, the fact 
that the zero momentum mode of the Goldstone propagator is identical to zero due to the rotation of 
the field ip as described in section has to be respected by excluding these modes in the summation. 
From the 1-loop propagators tt{p) and bb{p) we calculate the top and bottom time-slice correlator as 
defined in Eq. (j2ip via a Fourier transformation in the time direction. The physical fermion masses 
shown in Fig. [21d are then determined from the exponential decay of this correlator. 

From Fig. ^jp. one also learns that the finite volume effects are rather mild at A = 400 GeV with 
mH ■ Ls^t > 3.2 on the 32^-lattice, while the vev, and thus the associated cutoff A, as well as the 
Higgs boson mass itself vary strongly with increasing lattice size Lg at the higher cutoffs. In Fig. [3] 
we explicitly show this volume dependence of the vev and the Higgs boson mass, respectively. 

It is well known from investigations of pure 0(4) Higgs models [1^ |4^ [s^l that the vev as well 
as the mass receive strong contributions from the massless Goldstone modes, inducing finite volume 
effects starting at 0{Lj'^). The next non-trivial finite volume contribution was shown to be of order 
0{LJ^). In Fig. [3] we therefore plot the obtained data versus l/L^ and use the linear fit ansatz 

fiUL;')^A2l + Bi}l-L;' (32) 

to extrapolate to the infinite volume limit. The free fitting parameters ^i'm and -B^'L with the 
subscripts v and m refer to the renormalized vev Vr and the Higgs boson mass mn, respectively. 
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FIG. 3: (a) The dependence of the renormalized vev Vr — v/\^Zg on the squared inverse lattice side length 
l/Ll for the same model parameters used in the Monte-Carlo runs presented in Fig.[2K but for a greater set of 
different lattice sizes. The given hopping parameters k are associated to cutoffs ranging from A = 400 GeV to 
A ~ 1200 GeV. (b) The dependence of the Higgs correlator mass niH versus l/Ll for the same Monte- Carlo 
runs. The solid and dashed curves display the linear and parabohc fits as discussed in the main text. 



To respect the presence of higher order terms in 1 /L^ we include only the five largest lattice sizes 
in the linear fit. As a consistency check to test the dependence of these results on the choice of the fit 
ansatz, we furthermore consider the parabolic ansatz 

fv^mi-^s ) = ^iFm + ^v^m ' + ^v^rn ' (33) 

which we apply to the whole range of available lattice sizes. The deviation between the two fitting 
procedures are considered here as an additional systematic uncertainty for the infinite volume results 
determined from the linear fit in Eq. (|32p . The obtained fitting curves are displayed in Fig. [3^,b and 
the corresponding infinite volume extrapolations of the vev and the Higgs boson mass are listed in 
Tab.U 

So far, the presented results have been determined in the mass degenerate case, i.e. for yt — yb, 
which is easier to access numerically. This brings up the question of how the results are infiuenced 
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Linear fit 


Parabolic fit 


Final result 


K 




a(i) 


Aip) 


a(p) 


Vr 




0.12301 
0.12303 
0.12306 
0.12309 
0.12313 


0.1947(15) 
0.2761(24) 
0.3860(10) 
0.4883(03) 
0.6191(02) 


0.0535(20) 
0.0676(13) 
0.0826(15) 
0.0955(35) 
0.1063(11) 


0.1918(24) 
0.2766(21) 
0.3882(73) 
0.4903(03) 
0.6204(04) 


0.0517(15) 
0.0648(23) 
0.0825(14) 
0.0959(24) 
0.1090(17) 


0.1947(15)(29) 
0.2761(24)(05) 
0.3860(10) (22) 
0.4883(03) (20) 
0.6191(02)(13) 


0.0535(20)(18) 
0.0676(13)(28) 
0.0826(15)(01) 
0.0955(35)(04) 
0.1063(11)(27) 



TABLE I: The results of the infinite volume extrapolation of the Monte-Carlo data for the renormalized vev Vr 
and the Higgs boson mass mn- The extrapolation was performed by applying the linear fit ansatz in Eq. (|32|l 
and the parabolic ansatz in Eq. (|33} as described in the main text. The final results are taken from the linear 
fit. An additional, systematic uncertainty of the final results is specified in the second pair of brackets taken 
from the deviation between the linear and parabolic fit results. 



when pushing the top-bottom mass splitting to its physical value, i.e. mb/mt ~ 0.024. From Eq. (|22ll 
one expects the Higgs boson mass shift Smf^ to grow quadratically with decreasing yt and that is 
exactly what is observed in Fig. |4)d. Here, the bare top Yukawa coupling constant yt, the quartic 
coupling parameter, and the cutoff are held constant, while yi, is lowered to the physical ratio of yb/yt- 
The numerical data for the Higgs boson mass shift Sm^ are in excellent agreement with the predictions 
based on the effective potential, which are depicted by the dashed lines in Fig.HjD. However, the Higgs 
boson mass itself does not increase but decrease with decreasing yb as shown in Fig. [4^. This is because 
the first effect in the mass shift is over-compensated by the shift in the phase transition line, which 
is moved towards smaller bare Higgs boson masses m. This behaviour of the bare Higgs boson mass 
is demonstrated also in Fig. [Jb, where the numerical findings for the bare mass are compared to the 
corresponding expectations from the effective potential. Again very good agreement between both 
approaches is observed. 
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FIG. 4: (a) The Higgs boson mass versus the squared ratio of the top and bottom Yukawa coupling constants 
[yb/yt]^ on a 12^ x 32-lattice for constant cutoff A = 400 GeV, A = 0, Nf = 1, and yt fixed according to 
Eq. (|16|) . (b) The corresponding Higgs boson mass shifts Srnjj and bare masses versus [yt/yt]'^- The 
dashed curves represent the analytical predictions derived from the effective potential. 



Finally, we present the analytical infinite volume results, which are obtained from the effective 
potential in Eq. (|27p by replacing the sums with corresponding integrals. In Fig. [5] we compare these 
findings to the infinite volume extrapolation of the Monte-Carlo data taken from Tab. |T] and observe 
very good agreement. We also show the analytical infinite volume expectations for several other 
physical scenarios with varying values for Nf and yb/yt to study the influence of these parameters 
on the lower mass bound. The number of fcrmion generations (or here equivalently the number of 
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colours) as well as the bottom-top mass splitting have a strong impact on the lower mass bound. 
The solid curve in Fig. [5] with Nf = 3 and yb/ut — 0.024 comes closest to the actual situation in 
the Standard Model. However, direct Monte-Carlo simulations in this scenario on sufficiently large 
lattices would be too demanding with our available resources due to the huge mass splitting between 
the top and bottom quark. 




500 1000 1500 2000 2500 3000 3500 
Cutoff A [GeV] 



FIG. 5: The various curves show the cutoff dependence of the lower Higgs boson mass bound in the infinite 
volume limit derived from the effective potential for different physical setups {Nf,yi,/yt). These analytical 
findings are compared with the infinite volume extrapolation of the Monte-Carlo data taken from Tab. |T] 
depicted by the square symbols. The numerical results were obtained at {Nf — l,yb/yt = 1) and are in 
very good agreement with the corresponding dash-dotted curve. The solid curve represents the setup {Nf = 
'i^Vb/yt = 0.024) coming closest to the actual physical situation in the Standard Model. 



VI. OUTLOOK TOWARDS UPPER HIGGS BOSON MASS BOUNDS 

We now turn towards the determination of the upper Higgs boson mass bound m^(A). First, 
we investigate whether the largest Higgs boson masses are obtained at A = oo as expected from 
perturbation theory. This can indeed be observed in Fig. [6^ where we plot the physical Higgs boson 
mass m///a versus the quartic self-coupling constant in the strong coupling region of the model for 
a fixed cutoff A k, 1500 GeV. We compare these numerical findings for large but finite A with the 
corresponding A = oo result, which was obtained from our Monte-Carlo simulations by reparametrizing 
the Higgs field 4)x in terms of SU{2) elements. Since the integration of the molecular dynamics 
equations as needed in our Monte-Carlo algorithm becomes increasingly difficult with growing A, we 
present here only a few results on rather small 12^ x 32-lattices. However, within the available accuracy 
one can conclude from Fig.[6K that the largest Higgs boson masses are obtained at A = oo, as expected. 
We therefore derive the upper Higgs boson mass bounds in the following from simulations with infinite 
self-coupling constant. 

As discussed in the previous section the accessible energy scales A are restricted by Eq. p8)) . 
Employing a top mass of 175 GeV and a Higgs boson mass of below 700 GeV, as suggested by Fig. [5}:, 
we should be able to reach cutoffs between 1400 GeV and 2100 GeV on a 24^^ x 32-lattice. 

In Fig. we present the results for the top quark mass obtained in our Monte-Carlo calculations 
with degenerate bare Yukawa coupling constants fixed according to Eq. (flBl) and Nf = 1. Again 
we estimate the strength of the finite volume effects by rerunning the simulations on different lattice 
sizes. One sees that for the simulation runs with A < 2.1 TeV the volume dependence of the top quark 
mass has already reached its saturation plateau on the 24'^ x 32-lattices within the given accuracy. In 
comparison to the weak coupling region of the model the top quark mass depends significantly weaker 
on the cutoff. The discrepancy between the tree-level determination of the bare Yukawa coupling 
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constants and the targeted top quark mass is around 1% for the largest presented lattices and thus 
smaller than in the weak coupling region. However, in follow-up calculations the available data can 
be used for a fine tuning of the bare Yukawa coupling constants to reproduce the targeted fcrmion 
masses even more precisely. 

Finally, we present our results on the cutoff dependence of the upper Higgs boson mass bound in 
Fig. [6j:. Here, we decided to show only the propagator masses instead of the correlator masses. The 
reason is that the latter results are very noisy, since they are based on only those Lt Higgs modes 
in momentum space which have spatial momentum zero, and do - at the time being - not allow for 
a clear observation of the cutoff dependence of the Higgs boson mass. The propagator masses, on 
the other hand, are much more stable and do reveal the desired cutoff dependence of the upper mass 
bound. On the small 12^ x 32-lattices the Higgs boson mass seems to increase with the cutoff A. This, 
however, is only a finite volume effect. On the larger lattices the Higgs boson mass decreases with 
growing A as expected. Only when the cutoff becomes too high for the given lattice volume does the 
Higgs boson mass finally start increasing again due to the finite volume. Larger lattice sizes and more 
statistics are needed here to further follow the cutoff dependence of the upper mass bound. 
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FIG. 6: (a) The correlator Higgs boson mass versus the quartic coupling constant A in the strong coupling 
region of the model at A ~ 1500 GeV on a 12^ x 32-lattice. The dashed line depicts the corresponding A = 00 
result, (b) The cutoff dependence of the top correlator mass at A = 00. (c) The cutoff dependence of the Higgs 
propagator mass at A = 00. Finite volume effects are illustrated in panels b, c by rerunning the simulations 
on several lattice sizes. Runs with identical bare parameter settings are connected by lines to guide the eye. 
All presented results have been obtained in Monte-Carlo simulations with identical, degenerate bare Yukawa 
coupling constants fixed according to Eq. (|16|) and Nf — 1. 
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VII. SUMMARY AND OUTLOOK 

The eventual aim of our work is to establish non-perturbative upper and lower Higgs boson mass 
bounds based on first principle computations without relying on arguments such as vacuum instability 
or triviality. This can be achieved by investigating a lattice model of the pure Higgs- Yukawa sector 
of the Standard Model. The main idea of this approach is to apply direct Monte-Carlo simulations to 
determine the maximal interval of Higgs boson masses attainable within this model being consistent 
with phenomenology, i. e. with the physical values of the top and bottom quark masses and the vacuum 
expectation value of the Higgs mode. To maintain the chiral character of the Higgs-fermion coupling 
structure on the lattice we have considered here a chirally invariant lattice Higgs- Yukawa model based 
on the Neuberger Dirac operator. 

In the present paper we confirmed that the lightest Higgs boson masses are generated at vanishing 
quartic self-coupling, as expected from perturbation theory. For the mass degenerate case with equal 
top and bottom quark masses and A^^ = 1 we presented our numerical data for the Higgs boson 
mass obtained at A = for various cutoffs and lattice sizes. These results were compared to the 
analytical calculations of the Higgs boson mass based on the effective potential at one-loop order. 
We then extrapolated the numerical finite size results to infinite volume and compared them to the 
corresponding infinite volume effective potential predictions. In both cases very good agreement was 
observed. 

The main result of the presented findings is that a lower Higgs boson mass bound is a manifest 
property of the pure Higgs- Yukawa sector that evolves directly from the Higgs-fermion interaction for 
a given Yukawa coupling parameter. For growing cutoff this lower mass constraint rises monotonically 
with flattening slope as expected from perturbation theory. Moreover, the quantitative size of the 
lower bound is comparable to the magnitude of the perturbative results based on vacuum stability 
considerations [l9| . A direct quantitative comparison is, however, non-trivial due to the different 
regularization schemes accompanied with the strong cutoff dependence of the considered constraints. 

We also studied the Higgs boson mass dependence on the top-bottom mass splitting on some smaller 
lattices, where the calculations were feasible. Again the numerical findings matched the corresponding 
results based on the effective potential. The very good agreement between the numerical and analytical 
approaches justifies the presented effective potential calculations performed in the actual physical 
situation with Nf = 3 and mt/mi, w 40. A comparison with direct Monte-Carlo simulations in that 
setup on sufficiently large lattice sizes, however, seems to be too demanding for our available resources 
at the moment. We therefore have to rely on the analytical results in this case for the time being. 

Additionally, we turned towards the upper Higgs boson mass bound, which was obtained at infinite 
bare self-coupling after having checked that the heaviest Higgs boson masses are indeed generated at 
A = oo as expected. Instead of the correlator masses we decided to present the propagator masses at 
this point, since the latter results were stable enough to allow for the observation of the desired cutoff 
dependence of the upper Higgs boson mass bound in contrast to the correlator masses. Apart from 
finite volume effects the expected decrease of the upper bound with growing cutoff could clearly be 
seen. 

As a next step we will improve the statistics for the upper Higgs boson mass bound, extend the 
lattice volumes, and compare the numerical findings to corresponding calculations in renormalized 
perturbation theory, once sufficient precision is achieved. For the lower bound we plan to include 
higher order terms into the bare Lagrangian in order to investigate the universality of the obtained 
results. Taking these future improvements aside, our present results suggest a maximally allowed 
Higgs boson mass interval of 80Gev < tuh S TOOGev at a cutoff around Awl TeV. Eventually, we 
are also interested in studying the decay properties of the Higgs boson with this lattice model. 
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